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We present both numerical and analytical study of graphene roughness with a crys- 
tal structure including 500 x 500 atoms. The roughness can effectively result in a 
random gauge field and has important consequences for its electronic structure. Our 
results show that its height fluctuations in small scales have scaling behavior with a 
temperature dependent roughness exponent in the interval of 0.6 < x < 0.7. The cor- 
relation function of height fluctuations depends upon temperature with characteristic 
length scale of ~ 90 A (at room temperature). We show that the correlation function 
of the induced gauge field has a short-range nature with correlation length of about 
~ 2 — 3A. We also treat the problem analytically by using the Martin-Siggia-Rose 
method. The renormalization group flows did not yield any delocalized-localized 
transition arising from the graphene roughness. Our results are in good agreement 
with recent experimental observations. 
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I. INTRODUCTION 

An isolated two dimensional (2D) sheet of carbon atoms having a simple honeycomb 
structure is known as graphene^. According to the recent detailed measurements, these 2D 
zero-gap semiconductors reveal new features in their electronic properties. In particular, the 
low-energy quasi-particles of the system can formally be described by the massless Dirac-like 
fermions^. Doped graphene sheets are pseudo-chiral 2D Fermi liquids with abnormal effective 
electron-electron interaction physics 3 . From the experimental point of view, the melting 
temperature of thin films decreases with decreasing thickness and become unstable when 
their thickness reduces to a few atomic layers^. This observation supports the theoretical 
prediction that states: No strictly 2D crystal can be thermodynamically stable at finite 
temperatures-. Consequently, the expectation of not observing a free 2D material in nature 
lived on till the Geim's group discovered graphene^. Furthermore, because of the existing 
coupling between bending and stretching energy modes in any 2D material, one expects to 
observe a measurable undulations or at least very small roughness on the graphene sheet 
and will reduce electronic transport in graphene and affecting its thermal conductivity^. 

In this paper, we are interested in determining the temperature dependence of graphene 
roughness. Here roughness is defined as the variance of the height fluctuations in graphene 
due to temperature. Our numerical treatments suggest that its height fluctuations have a 
scaling behavior in small scale with a temperature dependent roughness exponent. We also 
determine the temperature dependence of the amplitudes of the height structure function. 
We show that the induced random gauge field has a short- range correlation function with 
correlation length ~ 2 — 3A. Finally, we formulate a field theoretical method to investigate 
the electron dynamics of the undoped graphene in such a random gauge field and show that 
there is no disorder- induced transition from delocalized to localized states 7 . This important 
result is in good agreement with the experimental observation that due to the graphene 
roughness localization is suppressed 8 . 

The rest of this paper is organized as follows. In Sec. II, we introduce an effective inter- 
action between carbon atoms that enters the molecular dynamic simulation to model the 
roughness of graphene . Section III contains our numerical calculations of graphene sim- 
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ulation and then using the latest results calculating the (3 function analytically by using 
the field theory of Martin-Siggia-Rose method. We finally conclude in Sec. IV with a brief 
summary. 

II. ROUGHNESS EXPONENT 

To study and arrive at a quantitative information of a graphene surface one may consider 
a surface with size L and define the mean height of the surface, h(L, A), and its roughness, 

L 

w(L, t, A), by the expressions like h(L, A) = ~ §\ h(x, A)dx and w(L, t, A) = (< (h — h) 2 > 
respectively. The symbol < ... > denotes an ensemble averaging. Here, A is an external 
factor which could be temperature in this problem - and t is the time. In the limit of large 
t, the roughness saturates and behaves as w(L, A) ~ L x<yX \ The roughness exponent, X-> 
characterizes the self-affine geometry of the surface. The common procedure for measuring 
the roughness exponents of a rough surface is to use a surface structure function, S(r) =< 
\h(x + r) — h(x)\ 2 > which depends on the length scale Ax = r. The surface structure 
function is equivalent to the statistics of height-height correlation function, C(r), and are 
related by S(r) = 2w 2 {l — C{r)) for stationary surfaces. The second order structure function 
S(r) scales with r as r^ 2 where % = £2/2— 

The atomic structure of graphene will force us to define two scaling exponents in x— (zig- 
zag) and y— (arm-chairs) directions. The exponents can be found via the second order struc- 
ture functions, S x (r x ) =< \h(x + r x ,y) — h(x,y)\ 2 > and S y {r y ) =< \h(x,y + r y ) — h(x,y)\ 2 >. 
The different scaling exponents in the zigzag and arm-chairs directions show the anisotropic 
nature of the roughness in graphene. We have used the empirical inter-atomic interaction 
potential, i.e., carbon-carbon interaction in graphite 12 , which has in addition three-body 
interaction for molecular dynamics simulation of graphene sheet to investigate its morphol- 
ogy and its dependence of roughness exponents on temperature. The two-body potential 
gives a description of the formation of a chemical bond between two atoms. Moreover, the 
three-body potential favors structures in which the angle between two bonds is made by the 
same atom. Many-body effects of electron system, in average, is considered in the Bren- 
ner potential, through the bond-order and furthermore, the potential depends on the local 
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environment. 

It is well known that the harmonic approximation resulting in bending instabilities due to 
soft wavelength phonons leads to crumpling of a membrane. It is important to note that the 
Brenner potential has anharmonic coupling between bending and stretching modes which 
prevents crumpling. We have considered a graphene sheet including a size of 500 x 500 
atoms with periodic boundary condition. Considering the canonical ensemble (NVT), we 
have employed Nose- Hoover thermostat to control temperature. Our simulation time step is 
1 fs in all cases and the thermostat's parameter is 5 fs. Therefore, we have found a stable 
2D graphene sheet in our simulation. 

III. NUMERICAL RESULTS 

In the top graph of Fig. 1, we have shown a snapshot of the graphene fluctuations at 
temperature 300A". Moreover, the magnetic field induced by the roughness is shown in 
the bottom graph to emphasize the randomness of the surface structure. Order of height 
fluctuations are about ~ 5A which is in good agreement with experimental observation^. In 
Fig. 2, we have plotted (in log-log scales), the structure function in the arm-chairs direction, 
namely S y (r y ), simulated at temperatures 10, 30, 100 and 300 K. Fig. 2 shows that for some 
characteristic length scales, the scaling behavior of the structure functions do not exist. The 
typical characteristic length scales are approximately 90A at room temperature which is in 
good agreement with experimental findings^ that is in the range of 50 — 100 A 

In Fig. 3, the temperature dependence of the scaling exponents, x' s ; i n both zigzag and 
arm-chairs directions are given. As shown at low temperatures, the exponents for zigzag- and 
armchairs- directions are about 0.7 (see also Ref.— ). However, the exponent for armchairs 
direction is greater than the one for zigzag direction at large enough temperatures. We 
have also used the bond order potential proposed by Ghiringhelli et al.—, and found good 
agreement between the results of the two potentials — . This is physically understandable, 
since our simulation has been performed for T < 700^ and the two potentials mainly differ 
at higher temperatures. Fig. 3 (inset) shows also the amplitude of the second moments, C x 
and C y , which are defined as, S x (r x ) = C x r x 2x and S y (r y ) = C y rl 2y in the scaling region, in 
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terms of temperature. 

To determine the characteristic length scales, we define the quantity Q(r) as the difference 
between the joint probability distribution function (PDF) of height fluctuations at two 
points. For instance, given y and y + r points, the Q{r) is calculated by P(h\,y\ h 2 ,y + r) 
and product of two PDFs, P(h%,y) and P(h 2 ,y + r]^-. Thus, 

Q(r) = J dh 1 dh 2 \P(h 1 ,y;h 2 ,y + r)-P(h 1 ,y)P(h 2 ,y + r)\ . (1) 

In Fig. 4 we have shown Q(r) as a function of r. This figure clearly indicates that the 
height fluctuations at scales of 90A and 125A for T = 300-ff and T = 30K, respectively, are 
almost independent. Obviously, the Q(r) becomes r— independent after these values of r 
(by considering its error bars). Furthermore, we have found the same value for correlation 
length for graphene with fixed boundary condition. As a consequence, this length scale is 
not an artifact of boundary conditions. 



IV. MARTIN-SIGGIA-ROSE EFFECTIVE ACTION 

Roughness of graphene results in a random gauge field and affects its electronic struc- 
ture^. The dependence of the hopping integral V on the deformation tensor is expressed 

by r = r + £-u lv where is given by u l3 = |{f| + fe + fefe + ftij*" Here ' 
Xi = (x,y) are coordinates in the plane and itj are the corresponding components of the 
displacement vector. In the presence of roughness, an effective Dirac Hamiltonian describes 
the electron states near the K-point, H = vpa(—iKV — -A), where vf = \/2>TQa/2h and 
A is the gauge field. The gauge field can be written in terms of the hopping integral V 
as: A x = ^(r 2 + T 3 - 2r : ) and A y = ^(T 3 - T 2 ) 19 . Labels 1,2 and 3 refer to 
the nearest neighbors atoms with vectors (— a/v3, 0), (a/2v^3, — a/2) and (a/2\/3, a/2), re- 
spectively. In Fig. 5, we have plotted the structure function of the gauge field A, namely 
5U( r ) — < |(A(x + r) — A(x)) • x| 2 > versus the scale in y-direction. The vectors r and x 
can be chosen in x— and y— directions. Therefore, we have four different types of structure 
functions. As shown in Fig. 5, the gauge field and the related magnetic field have small scale 
correlation with correlation length l c ~ 2 — 3A. We have checked the other three structure 
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functions of the induced gauge field and their cross correlation functions and find that they 
have small scale correlation and are almost statistically independent. 

In what follows we have developed a renormalization group analysis to investigate the 
Dirac equation in random gauge field and show that no delocalization-localization transition 
occurs for electrons in such a random gauge potential 20 ' 21 . 

The Lagrangian of Dirac fermions in (2+1) dimensions and in the presence of gauge 
potential is given by: C = i f dt f d 2 x-?/>7 M (<9 M — A^ip, in which 7 = a z ,rf = ic yi ^ 2 = 
—ia x , cr's are the Pauli matrices and 7's satisfy the Clifford algebra {7^,7^} = 2g fiu . The 
wave functions ip(x., t) and t/>(x, t) are 2D Dirac spinors and is a static random gauge 
field with a Gaussian distribution having zero mean value. The covariance is given by 
((Aj(x)A,-(x')) = 2D 5ij5{^x — x') where i,j = 1,2 and D is the intensity of its spatial 
fluctuations. This relation shows the spatially uncorrelated nature of the gauge filed. We 
note that the Fermi velocity of electrons in graphene is of the order of 10 6 m/s and the 
typical velocity of the height fluctuations is of the order of 30 m/s. Therefore, the random 
gauge potential will act as a quenched random field on electrons. Focusing on a single mode 
with energy e we get the following expression for the Lagrangian: 

C = j d^ii^dk - i-f k A k + e 7 )^ • (2) 
The expectation value of any operator O can be calculated as follows 

= fVrpDjOexpj-iC) ^ 
Z 

in which Z the partition function is defined by 

Z = J Vi/jV^O exp( J rf 2 x^(-7 fc 9 fe + 7 fc A fc - ie 7 °)V) 

= det(-7 fe «9 fc + 7 fc A fe -^7°) • (4) 

Introducing the Dirac bosons x an d X, one can re-expressing the above determinant as 
follows 

Z- 1 = J V X V X e W (- J d 2 *x{-l k d k + l k A k - *e 7 °)x) , (5) 
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and implicitly we have 



(O) = J VipVipVxVxOexpi J d 2 x^(- 7 fe ^+7%-ie 7 )V>- j d^xHl^k+^Ak-ie^x}. 

(6) 

Now by integrating the above result over the Gaussian variable whose probability 
density function is given by 



PocexpC-^i-y dxA|), (7) 

one reaches the following result for the averaging of expectation values over quenched random 
gauge 

«0» = J VipVi>VxT>xOex V (-S - S int ) , (8) 
where the free part of the effective action is 

So = J d 2 ^ k d k + iej )^ + j d 2 xX(l k d k + le^x , (9) 
and the interaction part 

S int = -D [ d 2 x [(^ + Xl k x)) ■ + Xl\)] , (10) 



where k = 1,2. The (3 function of the coupling D will determine its behavior under changing 
the scale. We found that in one-loop order the correction to the short-range roughness 
intensity (-Do) is proportional to k 2 , therefore this correction vanishes in long-wave length 
limit. One can see that the same will happen at higher order of of perturbation, leading 
us to the conclusion that the vanishing of the beta function in all orders should be the 
consequence of the Ward identity due the conservation of Dirac current which leads to an 
incompressible flow of electrons in the low energy limit ( e ~ 0). 

This result shows that the resistance against the electron flow due the interaction of Dirac 
fermions and roughness of the Graphene, remains unchanged under renormalization group 



flow towards large scales, and this in turn, excludes the possibi 



energy states (For more details, see for instance Ref[8.22.23 
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ity of the localization of low 



25|). 
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V. SUMMARY 



In Conclusion, we find the temperature dependence of the roughness exponents in dif- 
ferent directions of a graphene sheet by simulating the surface within molecular dynamics 
approach. We have used the Brenner empirical inter-atomic interactions for graphite which 
is a semiconductor. The point we are making in this publication is the further detailed ex- 
planations raised in Ref. 16j which roughness could affect. We answer the question whether 
the roughness could lead to localized electrons in graphene or not?. The correlation func- 
tion of height fluctuations shows that depending on the temperature, there are characteristic 
length scales in the order of ~ 90A at room temperature. We show that the induced gauge 
field has a short- range nature with correlation lengths approximately ~ 2 — 3 A. More 
importantly, roughness essentially can affect to the electronic properties like conductivity 
and modulation of the hopping integrals- 6 . We treat the problem analytically by using the 
Martin-Siggia-Rose method. The renormalization group flows do not yield any delocalized- 
localized transition due to roughness. In the present work, the effect of Dirac-like electrons 
on roughness are not considered. It would be of intrest to develop our work may use the 
quantum molecular dynamics simulation or ab initio Car-Parinnello molecular dynamics for 
Dirac-like electrons to investigate the dependence of Dirac-like electron on the graphene 
roughness. 
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FIG. 1: (color online) Snapshot of the graphene surface at 300K (upper graph). The sample 
contains a lattice size of 500 x 500 atoms. In the lower graph, we have plotted the induced 
magnetic field due to the roughness of graphene surface. 
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FIG. 2: (color online) Log-log plot of the second moment of height difference as a function of r, 
in y (arm-chairs) direction, which shows that for samples with temperatures 10, 30, 100 and 300 
K the hight fluctuations have scaling behavior in small scales. It indicates that the roughness 
exponent (the slope of the plots) decreases with temperature, and means that the surface will be 
rough at high temperatures. A similar figure can be found for the height-height structure function 
in x (zigzag) direction. 
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FIG. 3: (color online) The temperature dependence of the roughness exponents show that graphene 
is smoother in arm-chairs direction as compared to zigzag direction. For very smooth surface the 
exponent will be ~ 1.0. At high temperature limit the exponents approach a random noise exponent 
(i.e. 0.5). The anisotropy of the graphene is due to the fact that the lattice spacing in arm-chairs 
and zigzag directions are different. It shows that (inset) the amplitude of the second moments, C x 
and C y , are increasing function of temperature. 
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FIG. 4: (color online)Scale dependence of Q(r), defined by Eq.(l) as a function of r. It shows that 
the height fluctuations have characteristic scales of 90A and 125A for T = 300K and T = 30K, 
respectively. 
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FIG. 5: (color online) Log-log plot of the second moment of the induced gauge field (inset structure 
function of magnetic field) in y direction, for samples with temperatures 30, 100 and 300 K. The 
correlation length is about 2 — 3 A Similar figures can be found for the other three gauge field 
structure functions. 



